Estimate for the fractal dimension of the Apollonian gasket in d dimensions 
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We adapt a recent theory for the random close packing of polydisperse spheres in three dimensions 
[R. S. Farr and R. D. Groot, J. Chem. Phys. 131 244104 (2009)] in order to predict the Hausdorff 
dimension dA of the Apollonian gasket in dimensions 2 and above. Our approximate results agree 
with published values in 2 and 3 dimensions to within 0.05% and 0.6% respectively, and we provide 
predictions for dimensions 4 to 8. 

PACS numbers: 05.45.Df, 61.43.Gt. 61.43.Hv 



Leibniz [l[ first suggested that a plane area can be 
completely covered with discs, in an approximately self- 
similar manner, through a construction which involves 
starting with three equal touching discs, and then re- 
peatedly adding the largest possible disc which touches 
three neighbours, but does not overlap with any disc al- 
ready in the packing. The result is illustrated in figure 
[TJ According to Pappus of Alexandria, the problem of 
finding such osculating circles was first studied by Apol- 
lonius of Perga, in whose honour this 'Apollonian pack- 
ing' is named. A similar construction can be envisaged 
for spheres (where each added sphere touches four neigh- 
bours 2]). In higher dimensions, a construction based 
upon iterating the analogue of Soddy's formula or ap- 
plying iterated inversions [|| to hyperspheres will lead to 
overlaps @. Therefore in this paper we use the term 
'Apollonian packing' to refer to an 'osculatory packing' 
, which starts from d+1 equal, touching hyperspheres 
at the vertices of a regular d— simplex, and where repeat- 
edly, the largest possible d-dimensional hypersphere is 
added to the existing packing that does not overlap any 
already present. The added hypersphere touches d+1 
others at this stage (although this fact is not needed for 
the subsequent argument). 

Apollonian networks @ , which are graphs derived from 
Apollonian packings, have been suggested as models for 
real-world networks, such as social networks and hier- 
archical road systems 0]- In these contexts, networks 
based on Apollonian packings with any dimension, in- 
cluding d > 3 may be of practical relevance Q . 

In lower dimensions {d — 2,3), the physical signifi- 
cance of Apollonian packings is that they can be used as 
idealized models of high density granular materials, for 
example in high strength concrete Furthermore, re- 
lated constructions, such as space-filling bearings in two 
[Io| and three [Til ] dimensions, and random space filling 
bearings [l2| have been proposed as simplified models 
for the geometry of turbulence [13j or the broken mate- 
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FIG. 1: Apollonian packing of discs. 



rial near a geological fault. Random Apollonian packings 
of shapes other than hyperspheres have also been stud- 
ied [14|. In all these cases, the method of adding spheres 
is modified so that it is no longer the largest possible 
non-overlapping sphere which is added at each stage; for 
example in the case of bearings an additional constraint 
is needed to ensure a bichromatic colouring [111 ]. These 
modifications all have the effect of reducing the rapidity 
with which the packing approaches a volume fraction of 
unity as spheres are added and also alter the fractal di- 
mension of their residual sets. In recent work on random 
bearings [TBj the fractal dimension can even be varied 
continuously by choice of parameters. 

The residual set or 'Apollonian gasket' of such struc- 
tures is of practical relevance, since its surface area and 
volume (for the 3d case) are related to solvent adsorption 
and permeability to flow through the packing. These geo- 
metrical properties of the residual set are finite, provided 
the packing contains only spheres larger than a certain 
cutoff diameter. However, in the limit where spheres of 
arbitrarily small size are included, the residual set is frac- 
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tal in nature, and its Hausdorff dimension captures the 
essential geometric information [4|. 

Recent high-precision calculations have shown that in 
2d, the dimension of the Apollonian gasket is d,A,2 ~ 
1.30568 PJ], while in 3d, it is d A ,a « 2.4739465 0. ' 

Since the Apollonian packing is a special kind of sphere 
packing in three dimensions, then it is interesting to in- 
vestigate whether the recent approximate theory for the 
volume fraction of random close packings of polydisperse 
spheres, presented in Ref. [TtJ may shed some light on 
this problem also. The hope is that the essential geo- 
metric features of sphere packings captured in Ref. [17j ] 
might apply to non-random cases also. 

In the theory of Ref. [l7[ , we start with a known dis- 
tribution of sphere diameters P^d{D), where P^d{D)dD 
is the number fraction of the spheres with diameters in 
the range (D, D + d-D), and we ask what is the maximum 
random packing density which can be obtained? 

The procedure consists of several stages: First, P^d{D) 
is converted into a number distribution of one di- 
mensional rods P\d{L), by imagining a random non- 
overlapping (but not necessarily close packed) distribu- 
tion of spheres, passing a straight line through this dis- 
tribution, and counting each portion of the line within 
a sphere to be a rod. The resulting distribution of rod 
lengths is given [T3] by 



Pid(L) = 2L 



J™P 3d (D)dD 
/ °° P 3d (D)D*dD- 



(1) 



In order to simulate packing, we then imagine that this 
collection of rods interacts on a line through a hard pair 
potential which acts between each pair of rods Li and Lj 
through 



V(h) = 



oo if h < min(/Li, fLj) 
if h^minifLiJLj) 



(2) 



In Eq. ([2]) , h is the closest approach of the two ends of the 
rods, / > is a free parameter in the theory (which we 
explain later) , and the potential is able to reach through 
smaller rods which may be in the gap between the two 
rods under consideration. 

Finally, we search over all orderings of the rods, and 
find the ordering which occupies the maximum length 
fraction on the line. This search can be accomplished by a 
simple greedy algorithm [ItJ , and the final length fraction 
occupied by the rods is our estimate for the maximum 
random packing fraction of the spheres in 3d. 

As described, the theory depends on a free param- 
eter /, which in Ref. [17| is fixed by ensuring that 
the predicted close packing density for monodisperse 
spheres matches the known random close packing den- 
sity </>rcp ~ 0.6435. With this calibration, / sa 0.7654 
and the theory can be applied to arbitrary sphere size 
distributions. 

For the further development of this paper, we re- 
quire the generalization of this model to other dimen- 
sions. Therefore consider a polydisperse collection of 



1 



1/3 



1/9 



1/9 



FIG. 2: Analogue of Apollonian packing for rods on a line, 
using the potential of Eq. Q in the text with / = 1, and two 
iterations of fitting small rods into the gaps between large 
rods. 



d— dimensional hyperspheres, where P d d{D)dD is the 
number fraction of hyperspheres with diameters in the 
range (D,D + dD), and d > 2. If we consider passing 
a straight line at random through a single hypersphere 
of diameter D, then we will generate a collection of rods 
with a normalized length distribution given by 

Pu{L;D) = {d-l)LD 1 - d (D 2 -L 2 Y d ' 3 ^ 2 e(D^L), (3) 

where 6(x) is the Heaviside step function. 

The distribution of rod lengths generated from passing 
a line through a random distribution of d dimensional 
hyperspheres, will therefore be given by a convolution 
with P d d, but also taking into account that the collision 
cross section for the line with a hypersphere of diameter 
D scales as The result is 

/>oo 

P ld {x)oc / D d - 1 P dd {D)p ld {L-D)dD 1 (4) 

J D=L 

which with the correct normalization (obtained by re- 
versing the order of integration over D and L), gives the 
appropriate generalization of Eq. ([1} , namely 



Pid(L) = {d-l)L 



^{D 2 -L 2 )^y 2 P dd {D)dD 
J™Di-ip dd (D)dD 



(5) 



Now, consider the analogue of the Apollonian packing 
for rods on a line subject to the potential of Eq. (|2|). 
This consists of starting with a set of equal large rods, 
and then placing the largest possible rods into the gaps 
between them, which do not require the large rods to 
move. This process is then repeated iteratively, as in 
figure O 

In one unit cell of this structure, there is one rod of 
the longest length (which we take as unity) , which leaves 
a gap of size / to be filled by the smaller rods. In choos- 
ing and placing these smaller rods, we need 2° rods of 
length /(l + 2/) _1 , then in the remaining gaps, which 
are of length f 2 (l + 2 t f)~ 1 , we place 2 1 rods of length 
/ 2 (1 + 2/)~ 2 . Repeating this process, we have at itera- 
tion number j, V rods of length / 3+1 (l + 2f)-<< j+1 \ 

This implies that asymptotically, as L — > 0, the total 
number of rods of size greater than L behaves as 



/oo J 
P ld (L')dL'ocl + ^2\ 

i=0 

where L = f j+1 {l + 2f)-^ +1 \ Therefore 
P ld {L')dL' cx L x where x 



ln (iw) 



(6) 



(7) 
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Now, the distribution P\d{L) in Eq. (JT)) has a correspond- 
ing distribution P c u{D) of hyperspheres, which from Eq. 
§5§ is given asymptotically in the limit D — > by 



P dd {D')AD' oc D v where 



y 



1. 



(8) 



D 



We now link these results back to the dimension cLa of 
the Apollonian gasket in d dimensions in the following 
manner: According to Refs. [H, HH, the cardinal- 
ity of the set of spheres in an Apollonian packing, with 
curvature not exceeding k is given by 



N(k) oc K d 



(9) 



Combining Eqs. (J7]), (|8]) and ([9]), we therefore obtain 
our estimate for the Hausdorff dimension d A of the resid- 
ual set of the packing: 



d A 



-y = d-l- 



ln2 



In 



1+2/ 



(10) 



where / is the appropriate value for each dimension d of 
space. 

To complete the calculation, we need a value for the 
free parameter / in the theory. This will be done (as in 
the analysis of random close packing pj}) by calibrat- 
ing the theory for the monodisperse case. To obtain the 
rod distribution corresponding to monodisperse hyper- 
spheres, we use a collection of n ~ 50 000 rods sampled 
uniformly from the distribution of Eq. To do this, we 
take equal points on the inverse function of the integral 
of Eq. ([3]) , so that our rod lengths are given by 



Li = D 



1 

2n 



2/(d-\] 



-1 1/2 



(11) 



where i = 1 . . . n. 

Applying the greedy one-dimensional packing algo- 
rithm described in Ref. [17| , we calculate numerically the 
packing density of monodisperse hyperspheres as a func- 
tion of the parameter /, for each dimension d of space. 
The resulting curves for d — 2 to 6 are shown in figure [3] 

In order to apply the packing theory, we need the value 
of / for each dimension d. When applying the theory to 
random close packing in 3d, the calibration used was to 
ensure that the prediction for random close packing of 
monodisperse spheres was correct (l7j . 

For the Apollonian packing each hypersphere is added 
in such a way as to optimally fill the remaining avail- 
able space, and so the local geometry of packing will 
always be as efficient as possible. In order to capture 
this property, we choose / to give the maximum possi- 
ble local packing fraction of equal hyperspheres. This 
corresponds to placing equal osculating hyperspheres at 
the vertices of a regular d-simplex, and calculating the 
volume fraction occupied inside the simplex. We refer to 
this packing fraction as s imp , and it is illustrated for the 




FIG. 3: Plot of the maximum packing fraction pac k of d 
dimensional hyperspheres, as a function of the parameter /. 
The curves from bottom to top correspond to d = 2, 3, 4, 5 
and 6. 




FIG. 4: (a) The maximum packing fraction </> a im P for d — 2 is 
the ratio of the shaded area to that inside the complete regular 
2-simplex (equilateral triangle), (b) The same construction 
for d = 3. 



cases d = 2 and 3 in figure |U An alternative argument 
for this choice, is that we are calibrating / by using the 
true packing fraction of the first few hyperspheres in the 
Apollonian packing. In general, these will be of different 
sizes, but by taking the first d+ 1, we again only need to 
consider equal spheres at the vertices of a simplex. 

In 2d, this construction gives the same area frac- 
tion as a hexagonal packing, so </> S i mp = 71-/(2^3) « 
0.9069. In three dimensions, we find s i mp = 
3\/2 [cos" 1 (1/3) - tt/3] fa 0.7796, which is higher than 
can be achieved for any global packing of spheres in 3d 
(this limit is 4>{ cc = 7r/vT8 sa 0.74, achieved for the 
face centred cubic or hexagonal close packed arrangement 
20]). 

For higher dimensions, we calculate the simplex pack- 
ing fraction using a Monte-Carlo integration, noting that 
if one vertex of a regular simplex lies at the origin of 
d-dimensional Cartesian coordinates, then the other ver- 
tices can be chosen at the positions {s- 7 } where 

sj = (2d + 4 + AVTTdj [ (l + VTTd) Sij + 1 

(12) 

A point p chosen randomly (and uniformly) in (0, l) d can 
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TABLE I: Predictions for the Hausdorff dimension of the 
Apollonian gasket in d dimensions. The close packing density 
of spheres with centres at the vertices of a regular d-simplex is 
0sim P - The corresponding value of / from the packing theory 
is shown, along with the predicted Hausdorff dimension d A Te , 
and the actual Hausdorff dimension d Jf* if known 0, Gil • 



d 


0simp 


/ 


d A rcd 


7act 
"A 


2 


0.906900 


0.131025 


1.3060 


1.3057 


3 


0.779636 


0.394834 


2.4586 


2.4739 


4 


0.6478 


0.7864 


3.5848 




5 


0.5257 


1.325 


4.6840 




6 


0.4195 


2.047 


5.7603 




7 


0.330 


3.02 


6.8189 




8 


0.255 


4.35 


7.864 





be expanded as p = J2j <7j sJ > where qj = t^pi and 



{2d + A + Ay/TTd) 1 ' 2 

i + VTTrf 



The point p lies within the simplex if all the q^s are 
positive, and their sum does not exceed unity. We de- 
note the volume of this simplex by Ksi mp , which can thus 
be obtained by a Monte-Carlo integration, or from the 
analytic expression V s i mp = 2~ d / 2 y / (l + d)/d\. Further- 
more, consider the volume V sp h of that portion of a unit 
radius hypersphere lying within a large regular <i-simplex, 
when the hypersphere has its centre at one of the vertices 
of the simplex. The point p lies within this volume if all 
the q/s are positive, and J2j(lj) 2 < 1- Again, this allows 



us to calculate T4 p h stochastically. 

From these two quantities, the maximum packing frac- 
tion in a simplex is given by 



(bsimp = (d+l)V sph / (2 d V s[ 



(14) 



which is shown in table HI alongside the predicted val- 
ues of dA [from Eq. (fTtj|) ] and the actual values (where 
known) . 

From table HI we see that the predictions from this 
model in 2 and 3 dimensions agree with the known values 
to within 0.05% and 0.6% respectively, and predictions 
for higher values of d may be readily obtained. 

In conclusion, the packing theory of Ref. [13], which 
was designed to abstract the essential geometric features 
of random close packing, also appears to contain enough 
information to predict important features of the hierar- 
chical Apollonian packing. The extension of these argu- 
ments to more general Apollonian-type packings (such 
as space filling bearings [11] or random Apollonian pack- 
ings [141 ] ) will require further work, because the objects 
inserted into the packing are no longer maximal, which 
implies that both Eq. (J7J and the calibration of / will 
need to be modified. Nevertheless, we hope that further 
study of this or related theories will lead to more insights 
and further analytical results on both packings and gran- 
ular materials. 

Acknowledgments 

Fig. Q] was supplied by the user 'Time3000' on the 
'Wikimedia Commons' webproject. 



[1] G. W. Leibniz and B. Look, The Leibniz - Des Bosses 

correspondence Yale University Press (2007). 
[2] M. Borkovec, W. de Paris and R. Piekert, Fractals 2(4) 

521-526 (1994). 
[3] F. Soddy, Nature 139, 77-79 (1937). 
[4] B. B. Mandelbrot, The Fractal Geometry of Nature (W. 

H. Freeman & Co., New York, 1983). 
[5] H. Morall, J. Phys. A: Math. Gen. 27 7785-7791 (1994). 
[6] J. P. K. Doye and C. P. Massen, Phys. Rev. E 71(1) 

016128 (2005). 

[7] J. S. Andrade, H. J. Herrmann, R. F. S. Andrade and L. 
R. da Silva, Phys. Rev. Lett. 94(1) 018702 (2005). 

[8] Z. Z. Zhang, F. Cornelias, G. Fertin and L. Rong, J. of 
Phys. A 39(8) 1811-1818 (2006). 

[9] H. J. Hermann, R. M. Baram and M. Wackenhut, Brazil- 
ian J. Phys. 33(3) 591-593 (2003). 
[10] S. S. Manna and T. Vicsek, J. Stat. Phys, 64(3/4) 525- 
539 (1991). 

[11] R. Mahmoodi Baram, H. J. Herrmann and N. Rivier, 



Phys. Rev. Lett. 92(4) 044301 (2004). 
[12] R. Mahmoody Baram and H. J. Herrmann, Phys. Rev. 

Lett. 95, 224303 (2005). 
[13] G. Bachelor, Theory of homogeneous turbulence, Cam- 
bridge University Press, 1982. 
[14] G. W. Delaney, S. Hutzler and T. Aste, Phys. Rev. Lett. 

101, 120602 (2008). 
[15] P. G. Lind, R. M. Baram and H. J. Herrmann, Phys. 

Rev. E 77, 021304 (2008). 
[16] S. S. Manna and H. J. Hermann, J. Phys. A: Math. Gen. 

24, L481-L490 (1991). 
[17] R. S. Farr and R. D. Groot, J. Chem. Phys. 131 244104 

(2009). 

[18] D. W. Boyd, Mathematika 20, 170-174 (1973). 
[19] D. W. Boyd, Math. Comp. 39(159), 249-254 (1982). 
[20] T. C. Hales, Annals of Mathematics 162(3) 1065-1185 
(2005). 



